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1. Introduction 

Non-equilibrium quantum transport in fully-interacting systems is a barely-explored 
territory of theoretical physics. Equilibrium statistical mechanics of interacting systems, 
on the other hand, can sometimes be studied reliably by using powerful field-theoretical 
techniques (including conformal field theory and the Bethe ansatz). Bethe-ansatz integra- 
bility is useful in equilibrium even in the absence of scale invariance, permitting one to 
study the exact crossover behavior between critical points. However, integrability has for 
the most part been useful for calculating only thermodynamic quantities, excluding cor- 
relation functions and transport at non-zero temperature. However, by using the "quasi- 
particle approach" to integrability, we have recently shown that, quite unexpectedly, one 
can compute exact transport properties through a point contact, even out of equilibrium 
0. Here we discuss non-equilibrium transport through point contacts in a Luttinger liq- 
uid. This model is realized in resonant tunneling experiments through point contacts in 
quantum Hall effect devices 0. 

The key observation in [|I| was that (i) tunneling is integrable and that (ii) integrability 
defines a quasiparticle basis of the Hilbert space of the leads which is particularly suited 
to computing transport. In this basis, scattering processes at the point contact proceed 
without particle production. These quasiparticles are not free. However, their interactions 
can be incorporated exactly into non-fermi distribution functions, which govern the filling 
of single-particle levels (orbitals) with quasiparticles in the thermodynamic limit. Even 
though these distribution functions are not that of Fermi-Dirac, in an integrable model they 
can be computed exactly using thermodynamic Bethe ansatz (TBA) technology. Once this 
quasiparticle basis arising from integrability has been identified, non-equilibrium transport 
properties such as the exact I(V) curve and the conductance through the point contact can 
be computed exactly using a kinetic (Boltzmann) equation for these quasiparticles. This is 
possible even though the point-contact interaction is non-Gaussian, because the constraints 
of integrability give the exact S matrix of transmission and reflection amplitudes for these 
quasiparticles scattering off the point contact. This elastic single-quasiparticle S matrix 
is momentum-dependent, and, since there is no quasiparticle production in this basis, it 
is unitary. The S-matrix for scattering of a multi-particle state off the contact factorizes 
into a product of single particle S-matrices. This means that reflection and transmission 
processes of single quasiparticles by the point contact occur successively ("one-by-one"), 
and are the only scattering events. Therefore, these single quasiparticle processes describe 
the effect of the point-contact interaction fully and exactly. 
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The purpose of this paper is threefold: (i) We give the details of our exact non- 
equilibrium transport calculations through point contacts, using Bethe ansatz integrability. 
We have recently reported briefly on some of our results in [|J. (ii): We give universal 
explicit curves for the I(V) characteristic and the conductance of those point contact 
devices, for all temperatures, for comparison with future experiments, (iii): We prove the 
existence of an exact duality symmetry in this interacting system between weak and strong 
backscattering. In the context of the quantum Hall effect this corresponds to a duality 
between electron tunneling and Laughlin-quasiparticle tunneling. 

The Luttinger model is one of the simplest non-fermi-liquid metals It consists 
of left- and right-moving gapless excitations at the two fermi points in an interacting 
1-dimensional electron gas. In the past, this model had been difficult to realize exper- 
imentally. This is simply because in a one-dimensional conductor (such as a quasi-one- 
dimensional quantum wire so thin that the transverse modes are frozen out at low tem- 
perature), random impurities occur in the fabrication process. These impurities lead to 
localization due to backscattering processes between the excitations at the two fermi points. 
In other words, the random impurities generate a mass gap for the fermions. However, 
the edge excitations at the boundary of samples prepared in a fractional quantum Hall 
state should be extremely clean realizations of the Luttinger non- fermi liquids [ffl. These 
are stable because for 1/u an odd integer, the excitations only move in one direction on 
a given edge. Since the right and left edges are far apart from each other, backscattering 
processes due to random impurities in the bulk cannot localize those extended edge states. 
Moreover, the Luttinger interaction parameter is universally related to the filling fraction 
v of the quantum Hall state in the bulk sample by a topological argument based on the un- 
derlying Chern-Simons theory, and does therefore not renormalize. The edge states should 
therefore provide an extremely clean experimental realization of the Luttinger model. 

We study the tunneling conductance through a local backscattering potential in the 
Luttinger model, which gives a fingerprint of the non-fermi liquid state in the Luttinger 
metal. This situation is realized in resonant tunneling experiments through a point contact 
in v = 1 /3 quantum Hall devices . The point contact causes backscattering between 

right- and left-moving edge excitations, but since the coupling is only at a point, it can 
be controlled experimentally (and theoretically). The tunneling conductance can thus be 
viewed as a spectroscopy of the Luttinger non-fermi-liquid state in the quantum Hall edges. 
Indeed, the experimentally-measured linear-response conductance agrees remarkably well 
with our exact predictions based on the Luttinger model [IJ. This appears to provide 
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very convincing evidence that the Luttinger model describes the edge state in fractional 
quantum Hall devices. 

In this paper we study non-equilibrium transport through such an impurity in detail. 
Studying transport properties is crucial for making contact with experiment; thermody- 
namic properties such as the specific heat arising from a point contact are clearly not 
accessible experimentally. In the quantum Hall experiments one uses a 4-terminal geome- 
try of the quantum Hall bar, which is long in the x-direction and short in the y-direction. 
The left-moving (upper) edge of the Hall bar is connected to battery on the right such 
that the charge carriers are injected into the left-moving lead of the Hall bar with an equi- 
librium thermal distribution at chemical potential Similarly, the right- moving carriers 
(propagating in the lower edge) are injected from the left, with a thermal distribution at 
chemical potential /ir. The difference of chemical potentials of the injected charge carriers 
is the driving voltage V = (Xl — fiR. If V > 0, there are more carriers injected from the 
right than from the left, and a "source-drain" current flows from the right to the left, along 
the ^-direction of the Hall bar. In the absence of the point contact, the driving voltage 
places the right and left edges at different potentials (in the y-direction, perpendicular to 
the current flow), implying that the ratio of source-drain current to the driving voltage V 
is the Hall conductance. 

In the absence of any point-contact interaction, the source-drain current Iq(V) may 
be computed in a variety of ways (see e.g. [0,0). The resulting conductance is ve 2 jh 
(in linear response and at finite driving voltage V). When we include a point-contact 
interaction at finite driving voltage, more of the left-moving carriers injected from the right 
are backscattered than those injected from the left, resulting in a loss of charge carriers 
from the source-drain current. In this case we write the total source-drain current as 
I(V) = Io(V) + Ib(V), where Ib(V) is the (negative) backscattering current, quantifying 
the loss of current due to backscattering at the point contact. It is this backscattering 
current which we compute exactly here, for finite driving voltage V and all temperatures 
T. This computation is possible since backscattering does not deplete the (infinite) right 
and left reservoirs (the battery), so that we can use the individual thermal distribution 
functions for the left and right reservoirs. 

The I(V) characteristic as well as the conductance are described by three parameters, 
V, T and Tb, where Tb (analogous to the Kondo temperature in the Kondo effect) is 
a measure of the point contact interaction strength. However, due to the underlying 
quantum critical point, these observables are described by universal scaling functions of two 
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ratios, the scaled interaction strength Tg/T, and V/T. We display the exact differential 
conductance G(Tb/V,V/T) = 9I qX ' in fig. 1. As apparent from fig. 1, an interesting 
prediction of our exact solution is that the differential conductance exhibits a pronounced 
maximum as a function of the point contact interaction strength, whenever the ratio V/T 
exceeds a critical value (which depends on the filling u, i.e. on the Luttinger interaction 
constant). This is a pure non-equilibrium effect, since the maximum occurs only at finite 
driving V^O. Of course, the total current still decreases as the point-contact interaction 
strength is increased. This feature is consistent with the first-order perturbative results 
H for the I(V) curves at zero temperature, in the strong and weak backscattering limits. 

When there is a single relevant operator corresponding to the impurity, the model 
is integrable, so we compute the current and conductance exactly. There is only one 
relevant operator when v = 1/3, so integrability is generically observed without any 
fine tuning. Any other sample-specific details appear only in the irrelevant operators. 
Integrability allows the definition of a basis of massless charge-carrying 'quasiparticle' 
excitations @||. These quasiparticles are scattered one- by-one off the impurity with a 
momentum-dependent one-particle scattering matrix S of transmission and reflection am- 
plitudes. These amplitudes can be determined exactly by imposing the constraints of 
integrability, including the boundary Yang-Baxter equation and the boundary-crossing re- 
lation H. Furthermore, the quasiparticles are characterized by an thermal distribution 



function which can be calculated exactly using the thermodynamic Bethe ansatz [10,11 



This special behavior is a consequence of integrability, and it allows us to derive an exact 
rate (Boltzmann) equation for the conductance in this interacting theory. 

We should note that many of these results have been checked by using a completely 
different method of computation. Instead of the non-perturbative methods to be described 
below, one can study the model perturbatively in the interaction strength Tb- Closed- form 
results for all the perturbative coefficients of the free energy have been found; a simple 
technical assumption then also leads to closed-form expressions for all the perturbative 
coefficients of the conductance, even at finite V [fL2 |. There is complete qualitative and 
quantitative agreement between the two approaches. 

The outline of this paper is as follows: In sect. 2 we discuss how this system is realized 
in experiments on resonant tunneling through point contacts in fractional quantum Hall 
devices. In sect. 3, we map the problem into two decoupled theories, one of which is 
affected by the backscattering interaction, and another one which is not. In sect. 4 we find 
the exact quasiparticle spectrum, exact S matrices for quasiparticles scattering among 
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themselves and off of the boundary, and the exact thermal distribution functions, using 
the fact that the interacting theory is integrable. In sect. 5, we derive an exact (Boltzmann) 
equation for the current in terms of the S matrix and the distribution functions. This gives 
us equations for the exact conductance. These can easily be solved numerically, and we 
present curves for a variety of values. At zero temperature, the equations simplify and we 
present more explicit analytic results in sect. 6. This enables us to derive an exact duality 
between strong and weak backscattering limits, out of equilibrium. 

2. The fractional quantum Hall effect 

2.1. Experimental setting for resonant tunneling 

Experiments on resonant tunneling between two v = 1/3 edges have recently been 
performed by Milliken, Umbach and Webb [0. We briefly review the experimental setup 
schematically (for details see A fractional quantum Hall state with filling fraction 

v = 1/3 is prepared in the bulk of a quantum Hall bar (discussed in Section 1). This 
means that the bulk quantum Hall state is prepared in a Hall insulator state (longitudinal 
conductivity o~ xx = 0), and that the (bulk) Hall resistivity is on the v = 1/3 plateau where 
o~ xy = (l/3)e 2 //i. This is achieved by adjusting the applied magnetic field, perpendicular 
to the plane of the bar. Since the plateau is broad, the applied magnetic field can be varied 
over a significant range without affecting the filling of v = 1/3. Next, a gate voltage V g is 
applied perpendicular to the long side of the bar, i.e. in the y direction (see Section 1) at 
x = 0. This has the effect of bringing the right and left moving edges close to each other 
near x = 0, forming a point contact. Away from the contact there is no backscattering (i.e. 
no tunneling of charge carriers) because the edges are widely separated, but now charge 
carriers can hop from one edge to the other at the point contact. 

The linear-response source-drain conductance as a function of temperature and gate 
voltage V g has been measured experimentally @ . As the gate voltage is swept through, the 
conductance signal shows a number of resonance peaks, which sharpen as the temperature 
is lowered. These resonance peaks occur for particular values V g = V* of the gate voltage, 
due to tunneling through localized states in the vicinity of the point contact. Ideally, on 
resonance, the source-drain conductance is equal to the Hall conductance without point 
contact, i.e. G resonance = (l/3)e 2 /h. This value is independent of temperature, on reso- 
nance. Now, measuring the linear response conductance as a function of the gate voltage 
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near the resonance, i.e. as a function of 5V g = V g — V*, at a number of different temper- 
atures T, one gets resonance curves, one for each temperature. These peak at SV g = 0. 
Scaling arguments imply |J that those experimental conductance curves should collapse 
onto a single universal curve when plotted as a function of 5V g /T 2 / 3 . Indeed they do col- 
lapse quite well 0. Moreover, the resulting universal curve should be a unique signature 
and fingerprint of the v = 1/3 edge state in the leads connected to the point contact. 

2.2. Comparison of the linear response resonance lineshape with our exact Luttinger 
model predictions 

To make contact with the Luttinger model, we state (as will be explained below 
in more detail) that the experimental parameter SV g should be related to the Luttinger 
backscattering interaction by 5V g oc T^ 3 . In particular, at the resonance value Tb = 0, 
there is no backscattering at all. 

We have compared [|1| our exact predictions for the linear response conductance scaling 
curve with the experimental data [0 as well as Monte Carlo calculations in 0. The 
agreement between the Monte Carlo simulation and our exact scaling curve is excellent. 
The exact value of the universal parameter K (defined so that G(X) = KX~ 6 for X large 
and G(X) = (1 - X 2 )/3 for X small) is K = 3.3546... (where X « .74313(T B /T) 2 / 3 ). 
(The value K « 2.6 quoted in seems to have been slightly underestimated there.) 

The comparison with the existing experiments by Milliken et al is not completely 
straightforward, since the conductance at the resonance peak in the experimental data 
decreases with temperature and is well below its resonance value e 2 /3/i. This difficulty 
arises since in order to achieve the resonance condition in the Luttinger model, two param- 
eters need to be tuned, since the point contact will in general not possess reflection-parity 



symmetry about the point x = In the experiments performed so far only one 

parameter, namely V g , has been tuned. For that reason, the conductance peaks do not 
have their maximum height G resonance = (l/3)e 2 //i, but are smaller, and, furthermore, 
the peak height does decrease with temperature, reflecting the fact that the experimental 
peak is not a perfect resonance. This problem can be remedied in a future experiment, 
by varying two parameters, namely V g and the magnetic field on the plateau, instead of 
only one parameter, to achieve resonance. Nevertheless, even when only the gate voltage 
is tuned to resonance, the experimental data for the conductance signal as a function of 
the gate voltage and temperature collapse well onto single scaling curve. Thus it makes 
sense to compare the this experimental curve with our exact conductance curve, computed 
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from the Luttinger model. The agreement is quite good, given the large scatter of the 
data in the tail of the resonance curve. In particular, the data clearly show the predicted 
G oc T 4 /(V g - V*f behavior in the tail. 

2.3. Predictions for future non- equilibrium transport measurements 

In principle, there is no reason why the above mentioned measurements could not 
be extended to finite driving voltage V [one should not confuse the driving voltage V, 
which is the difference between the chemical potentials between the injected left- and 
right-moving charge carriers, with the gate voltage V g , which gives rise to the coupling 
constant Tb in the Luttinger liquid theory]. So far, only the linear-response conductance, 
G = limy^o I(V)/V has been measured. More generally, one could attempt to measure 
the conductance at finite driving voltage V. Note that our exact results predict the shape 
of the universal scaling function 

G{T B /T,V/T), T B = C(5V g f/ 2 

as a function of two ratios. The non-universal parameter C is determined by fitting the 
experimental data to the universal curve; this is the only unknown quantity. The linear- 
response conductance is the limit of this function as V — > 0, where it becomes a function 
only of one ratio Tg/T. Note that the conductance at finite driving voltage V also describes 
a resonance lineshape. A particularly interesting feature of the conductance at finite driving 
voltage V is seen in fig. 1. The conductance as a function of T B oc (SV g ) 3 ^ 2 develops a 
pronounced peak when the ratio of driving voltage to temperature eV/ksT exceeds a 
critical value 7.188.... For a typical low temperature T = 50mK used in the data of 0, 
this would correspond (10k bK = ImeV, or, = 1) to a driving voltage of V* ~ 35// V. 

The current data were taken at an 'excitation voltage' of 1/tV, which corresponds to a ratio 
eV/ksT ~ 0.2. We have plotted the exact results for the aforementioned scaling functions 
in Figs. 1 and 2. One sees clearly from these plots that eV/ksT ~ 0.2 corresponds to the 
linear response regime. Perhaps the most significant feature displayed in Fig. 2 is a very 
dramatic non- equilibrium broadening of the resonance curve for values of the ratio V/T 
even well below the occurance of the maximum. In terms of numbers, the curve broadens 
by a dramatic amount already at eV/ksT ~ 2 or 3, well before the onset of the maximum 
(at eV/ksT « 7.188). This broadening should be easily visible experimentally, since it 
would correspond to an excitation voltage V ~ 15fiV, just a factor of 15 higher than the 
ones used by Milliken et al., at T = 50mK. 
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Clearly, at larger voltages one will have to worry about larger currents flowing in 
the sample, which might lead to possible complications. These complications would for 
example include non-universal effects arising from sample heating which are clearly not 
included in the Luttinger theory and which could mask the underlying universal Luttinger 
non-fermi liquid physics. However one should notice that for most part of the conductance 
curve the conductance is very small, implying flow of very small currents, so that sample 
heating should not appear to be a problem. 

Futhermore, it is important to notice that the scaling function depends only on the 
ratio V/T. This means the same critical ratio (V/T)* = 7.188 could also be achieved by 
working at lower temperatures. It might be reasonable experimentally to perform mea- 
surements at lower temperatures. This would reduce even more the values of the driving 
voltage where the onset of the universal non-equilibrium features (the most dramatic one 
being the broadening) become clearly visible. 



3. Bosonizing and mapping to an integrable model 

The Luttinger model is the most general model of a single massless fermion in one 
dimension. The fermion interaction is governed by a single parameter gLutt, so that when 
gLutt > the interaction is repulsive and gLutt < it is attractive. In the absence of the 
impurity, the Hamiltonian is 

H = n J dx{j 2 L + J 2 R + g L uttJLJ R ) (3.1) 

where the left movers are coupled to the right. We use the well-known map of this model 
to free massless left and right moving bosons $l and | H}| . To do so one bosonizes the 
currents Jl = — j^d x $>L and Jr = -^d x $R. The fermion operators take also the form 
^11 =: e l ® R : and ^l ='■ e l ® L :. By introducing new fields [15] using the transformation 
(canonical up to a global factor y/a) 

a + 1 a — 1 

^ a-1 a + 1 (3 ' 2) 
®l =^—¥r h — y~ </?i ' 

one can decouple the interaction in the bosonized hamiltonian that reads then 

h = 1 J dx [jl + f R ] = JL J n 2 + (a^) 2 , (3.3) 
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where IT is canonical momentum conjugate to <f>. The coefficient v is related to the old 
Luttinger coupling via g Lutt = -2g^j, \ = jf^i- The currents j L = -^d x (p L and 
3r = ^d x (p R are then normalized so that (^(^1)^(^2)) = (27r )2(^ 1 _ X2 )2 which makes 
the meaning of v as a chiral U(l) anomaly explicit. The operators ipR =: e^ i ' v ^ <pR : and 
ipL =: e( l t u ^ L : are sometimes referred to as "Luttinger hyperfermion" operators when 
1/u is an odd integer |16|]. There are two U(l) symmetries of this model corresponding to 
left and right (hyper) fermion number: we introduce the corresponding conserved charges 
Ql = f, dxjL, and Qr = J, oIxjr. The corresponding "hyperfermions" have charge e as 
do the original electrons, as can be seen by computing their commutators with the charge 
operators. 

We now include an impurity, which we assume is localized at the origin x = 0. It 
couples left- and right-moving fermions with the interaction ^ l (0)^r(0) + \I/^(0)\I/ l(0). 
Such a coupling destroys the separate conservation of the two charges Ql and Qr: only the 
total charge Ql + Qr remains conserved, corresponding to simultaneous phase translations 
of (fL and (fR. In a renormalized effective theory, all terms that are not forbidden by total 
charge conservation symmetry can appear in the hamiltonian. These may be represented 
in terms of the bosons <£>l, (Pr by a backscattering contribution to the Hamiltonian MM : 



H 



B 



1 ^^ e inipL(x=0) e -inip R (x=0) _|_ e -in<p L (x=0) e inip R (x=0) j ^3 ^ 



When we describe the Luttinger theory in terms of the bosons (pi, and (fR, i.e. after the 
Luttinger bulk interaction gLutt has been disentangled, this term describes hopping of 
quasiparticles of charge ve (the Laughlin quasiparticles in the case of FQHE) between the 
left- and right-moving edges. Only the terms with n 2 < 1/v are relevant. For v > 1 all 
terms are irrelevant; this case is more appropriately described by a "dual" picture [|J. At 
v = 1 there is a single exactly marginal operator. This preserves the conformal invariance 
of the fixed point, so conformal techniques are applicable, making possible the calculation 
of the complete partition function [jlTj. We now focus on the case v < 1. For the model 
to be integrable, it seems that only one relevant perturbation is allowed 0. Thus when 
1 > v > 1/4, the backscattering interaction is automatically integrable without any fine- 
tuning. The experimentally- measured value of v = 1/3 falls into this range of parameters. 
For v < 1/4, the model has to be fine-tuned in order to be integrable. In particular, for 
1/9 < i/ < 1/4 a single parameter (A2) needs to be tuned to zero in order to achieve 
integrability. In general, we will be treating the model with A n = for n > 2. 
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In order to transport net charge through the impurity, we place the injected left 
movers and the right movers at different chemical potentials. This amounts to adding a 
term e{Qi — Qr)V/2 to the Hamiltonian describing the thermal weighting of the injected 
charge carriers. Notice that even though this problem is out of equilibrium for non-zero 
V, the charge carriers injected into the leads are thermally weighted with an equilibrium 
distribution function corresponding to or fiR. Thus we are studying the coupling 
between two equilibrium distributions. 

It is convenient to introduce the "backscattering temperature" Tg, which is the scale 
generated by the relevant point contact coupling constant at one point: 

T B = CtX 1 ^ 1 -^ 

where C\ is some non-universal constant. When there is no backscattering, Tb = 0. The 
exponent follows from a simple perturbative analysis 0. For v < 1/2 and T > there 
are actually no divergences in the coefficients of the perturbative expansion; we exploit 
this fact in Sect. 6 to derive C\ exactly. The conductance G and current / are universal 
scaling functions of two dimensionless ratios, for example Tg/T and V/T. Since Tb is 
a scale introduced by the (local) impurity, a bulk quantity like the leading part of the 
thermal distribution functions is independent of the point contact interaction and can be 
computed without the point contact. In other words, the corrections vanish with the size 
of the system. 

We now map the model including the impurity to an integrable model. By taking 
(non-local) linear combinations of the bosons, we can first map the left- and right-moving 
bosons to a purely left-moving system. Calling these left movers even and odd, we have 

f (3-5) 
4>°{x + t) =-j=[ip L {x : t) -<p R (-x,t)] 

where these particles are defined on — I < x < I. One replaces jliJr with j e ,j° in the 
Hamiltonian ( |3.3|) , where j e ^°(x + t) = (l/y/2)\ji J (x,t)±jR(—x,t)] are the charge densities 
of even and odd bosons. For purposes of computing the backscattering current, the non- 
locality of this transformation does not matter (it clearly would if we were to compute 
spatially-dependent Greens functions). Two major simplifications arise from this change 
of basis. First, the interaction involves only the combination <^l(x = 0) — <pr(x = 0) so 
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only the odd boson (jf interacts, while the even boson <p e remains free [ffBJ. (This was 
the crucial step in solving the X-ray edge problem in a Luttinger liquid Second, the 

backscattering current Ib can be entirely expressed in terms of the odd boson theory. This 
is easy to see: The even and odd charges are related to the charges of the original left- and 
right-moving edges by AQ = Ql — Qr = \^Q° and Ql + Qr = V%Q e - Thus Q e measures 
the total charge on both edges and is conserved even in the presence of the interaction. 
Therefore, the even boson decouples from the quantities we study (although it does enter 
into the determination of the spectrum at the fixed points jn| ) and will play no role in the 
following. Moreover, the backscattering current is the rate at which the charge difference 
between right and left edges decreases (see below). It is thus directly related to the odd 
charge by dtAQ. 

The constraint of integrability is a very powerful one. As we will see, it enables the cal- 
culation of exact thermal distribution functions in an interacting theory. Not surprisingly, 
only certain models are integrable. The impurity is also significant; only particular types 
of impurity couplings preserve the integrability. The Luttinger model without the impurity 
is well known to be integrable, but it is easily shown that before the non-local map ( |3.5|) , 
the impurity destroys the integrability. However it is integrable in the odd-boson basis. 
The integrability of the impurity interaction in the Luttinger model was first established 
by || in the context of the boundary sine-Gordon model. In order to make contact with 
this work, we map our theory involving the massless left-moving odd boson (jf on the 
line, with the point contact interaction at the origin x = 0, into a boundary sine-Gordon 
problem. This is done using a standard "folding" procedure. From the left-moving odd 
boson <p°(x + 1) on the full line, we define right and left moving odd bosons on the half 
line x > by 

4>° L {x,t) = (f)°(x + t), x>0 

4>° R (x,t) = (j>°(-x + t), x>0 

This model is no longer chiral, but lives on the half-line. 

There are many field theories on the half-line which are integrable, the Kondo model 
being a celebrated example ||20|| . Recently much effort has gone into understanding their 
properties. In the following we exploit the results of ||. The odd-boson model with 
interaction becomes the massless limit A —>■ of the massive sine- Gordon model on the 
half-line 

H SG = J- [ dx \(d x <j) ) 2 + (IT) 2 + Acos y/2<j>°\ + Ai cos -^°(0) (3.6) 
Sins Jo i J V2 
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[11° is the canonical momentum conjugate to the odd-boson field This model is inte- 
grable for any value of the bulk mass A, and the Luttinger parameter v [||. However, the 
factor of two ratio in the argument of the bulk and boundary cosines seems to be necessary 
for the integrability. In the standard sine-Gordon conventions where the first term has a 
1/2 in front, our normalization corresponds to /3§ G = 8ttv. A lattice regularization of this 
field theory, the XXZ spin chain with a a x perturbation on the boundary, was long ago 



shown to be integrable [21 



Notice that we can "fermionize" to get back a Luttinger model (for A^O we get the 
massive Thirring model) on the half-line for the odd-boson degree of freedom. However, 
because of the non-linear change of basis in ( |3.5| ), the properties of the odd-boson Luttinger 
model are not the same as the original model. In particular, the original Luttinger model 
is a free fermion when v = 1, while the odd-Luttinger fermion is free when v = 1/2. This 
shift arises because of the \pi in (|3.5| ), which is necessary to keep the same normalization 
of the kinetic term. Thus an interacting fermion on the full line can be mapped to a 
non-interacting one on the half-line. This provides a simple way of understanding the 
results of J22|, where the v = 1/2 model is mapped onto two Ising models, one of which 
has a boundary magnetic field. Two Ising models are well known to be equivalent to a free 
Dirac fermion (up to boundary conditions), and it is easily seen that Hb corresponds to a 
boundary magnetic field on one of them. 

To recover our massless model, we simply have to take the A — > limit @. So we see 
that, after the few mappings described above, our impurity problem is integrable as long 
as only one of the impurity coupling constants X n in ( |3.4|) is non-zero. As noted before, 
this is natural for 1 > v > 1/4, where power counting shows that only the first coupling 
constant Ai is relevant. 



4. The quasiparticles, their S matrix and their distribution functions 

4-1. The quasiparticle spectrum 

We find the appropriate quasiparticle basis by studying the model (|3lJ) at arbitrary A. 
As discussed carefully in |J, the quasiparticle spectrum remains the same in the massless 
limit A — > of the massive sine-Gordon model on the half-line ( |3.6|) . Introducing the 
parameter A and then setting it to zero is not necessary for solving the problem, but it 
gives an easy way of finding the basis where the quasiparticles scatter off the boundary 
without particle production. 
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Integrability means that there is an infinite number of conserved quantities which 
commute with each other and with the Hamiltonian, even in the presence of the impurity 
(EOD B • From this we see already that in an integrable system a particular basis of Hilbert 
space, in which all the infinite conserved quantities as well as the hamiltonian are diagonal, 
plays a special and simplifying role. It is this basis on which we will focus. This basis has 
a "quasiparticle" structure similar to the "particle" Fock space of a non-interacting theory. 
In particular, the eigenvalues of the infinite conserved quantities in this basis have the form 
Y2i Pii where pi are momenta of individual quasiparticles, and where n runs over an infinite 
subset of the positive integers (each n labels one conservation law). In the Luttinger model 
these run over n odd. The identification of pi as momenta of "quasiparticles" arises from 
the form of the eigenvalue of the Hamiltonian in this basis, which is H = vf ^2iPi- This 
means that the energies of the "quasiparticles" are additive, justifying the association of 
the particle concept with those eigenstates. 

These conservation laws have important consequences for scattering. The quasiparti- 
cles of this basis must scatter off the impurity without particle production, i.e. one-by-one. 
This means that the scattering matrix off the point contact is a product of 1-body S matri- 
ces, one for each quasiparticle. Away from the impurity (in the "bulk"), the quasiparticles 
scatter off of each other with a completely elastic and factorizable two-body scattering ma- 
trix S bulk . This follows from an old kinematic argument |f23| , |24|| which applies when there 
is any conservation law with n > 1. Factorizability means in the bulk that the iV-body 
bulk S matrix is a product of 2-body S matrices S . Completely elastic means that 
individual momenta are conserved in a collision: all that can happen is that the momenta 
of the particles get permuted. This does not mean that bulk scattering is trivial; internal 
quantum numbers can change, and even if scattering is diagonal (i.e. internal quantum 
numbers do not change), the particles can have a phase delay, i.e. the S matrix may be a 
momentum-dependent phase. 

There are many bases for the Hilbert space of the massless odd-boson theory, which 
are related by not-necessarily-local mappings. For example, plane waves obviously are 
eigenstates of Hq, but they are not eigenstates of Hq + Hb ||. The basis of particles 
of ( |3.6| ) which are eigenstates has been known for some time [ |24| . (Indeed, in a massive 
theory there is only one particle basis; it is only for A = in ( |3.6| ) that there is a choice.) 
At any value of the spectrum contains a kink (+) and an antikink (— ). These carry 
(odd) charges Q° = 1/V2 and —l/y/2, respectively. Moreover, for n — 1 < 1/u < n, there 
are n — 2 "breather" states which have no charge. These breathers exist in the regime 
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where the fermion interaction in the odd-Luttinger model is attractive; in this language 
they correspond to fermion-antifermion bound states. These particles span the Hilbert 
space of the left-moving odd boson; we label them by indices j, k, . . . running over the 
kink (+), antikink (— ) and breathers (b). One can in fact check explicitly that that these 
particles are the solutions of the classical limit of (|3.6|) for any value of A 0. 



Henceforth, we set A = in (|3.6|) so that the particles are massless. We also ignore 
the even boson, since it does not affect the current or conductance. We will find more 
convenient to use the unfolded language, so we will continue to discuss a purely left-moving 
theory on the full line with an impurity. 

Since the particles are massless, a left mover has dispersion relation E = —p. Instead 
of momentum, we use rapidity 9, which for a particle of type j is defined as 

E = —p = rrije 9 , 

where irtk = M sin(/c7r^/2(l — u)) for the fcth breather and m± = M/2. The overall scale 
M is arbitrary and cancels out of all physical results. 

Notice that the momenta of the quasiparticles all have one sign (the rapidity 9 is real) . 
Since this might be unfamiliar, we now express the charged fermions occurring for v = 1/2 
in terms of kink and antikink quasiparticles; for v = 1/2 there is no breather. Recall that 
for v = 1/2 we can re-fermionize the odd-boson theory. Thus we obtain a single non- 
interacting, spinless charge fermion tp^ (p),ip(p) where the momentum p runs over all real 
values (positive and negative). This fermion satisfies canonical anticommutation relations, 
since it is non-interacting. The kinks and antikinks for this simple theory can be defined 
using a canonical particle-hole transformation: 

tp+ip) = ip{p) (kink) 
ip-{p) = V^( — p) (antikink) 

for p > 0. The left-hand-side defines kink and antikink annihilation operators. Similarly, 
the kink and antikink creation operators are the hermitian conjugates: 

ijj\_(p) = if)\p) (kink) 
ip_(p) = 4>(—p) (antikink) 

for p > 0. Thus for v = 1/2 kinks and antikinks are just particle- hole transforms of 
ordinary fermions, and have only one sign of momentum. Unfortunately, for v ^ 1/2 the 
mapping is not as straightforward. 
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While there are no statistics in 1 + 1 dimensions (there is no way to bring a particle 
around another without interacting), a crucial issue is how many particles are are allowed 
to occupy each level. The answer is well established in the massive sine-Gordon model: 
only one quasiparticle is allowed per level (they are like fermions). This is basically because 
the Bethe-ansatz wave function used to solve the model or its regularized lattice versions 
vanishes when two or more excitations have the same rapidity. It is possible to study 
similarly a regularized lattice version of the massless limit of the sine-Gordon model [ p5| , p6| 
and to see that the same exclusion principle holds in that limit. This does not contradict 
other ways of describing the excitations, where the particles are free but obey "exclusion 



statistics" [p7j . These particles seem to be related to the quasiparticles of the Calogero- 
Sutherland model, or the Haldane-Shastry spin chain. In the realization of the Luttinger 
liquid in the fractional quantum Hall effect discussed in sect. 2, these particles are the 
Laughlin quasiparticles restricted to the edge. These particles are clearly not the same as 
ours; for example they are the quasiparticles of the original Luttinger impurity problem, 
while ours are those of the odd boson theory, which is found from the original by a non- 
local map. This means that there are (at least) two sets of bases of quasiparticles for the 
Luttinger liquid, one where the particles are free but fill levels in a peculiar manner, and 
another, where the particles fill levels like fermions but interact. In the first approach, 
however, it is not known yet how to include an impurity. If this problem is solved, it would 
be very interesting to compare the two methods of calculation. 

4-2. The bulk and impurity S matrices 

These massless kinks and breathers interact with each other and with the impurity. 
Even though the model away from the impurity is a free boson, an effect of using this quasi- 
particle basis is that the particles interact even in the bulk. We describe these interactions 
with bulk and impurity scattering matrices. First, we have an S matrix for a particle 
scattering off of the impurity. Since the integrability requires that the particles scatter 
one-by-one, this can be described by one-particle S matrix elements. Kinematically, all 
that happens is the particle goes through the impurity with a phase delay. Because the im- 
purity interaction violates odd-charge conjugation, it is possible for the positively-charged 
kink to scatter into a negatively- charged antikink when going through the impurity. 

Dimensional analysis tells us that the S matrix elements must depend only on the 
ratio p/Ts- If we define the "backscattering rapidity" by the relation 



M .. 
T B = —e° B . 
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this means that all impurity S matrix elements depend on the rapidity difference 9 — 9b- 
Thus the impurity-kink S matrix consists of the elements S ++ (9 — 9b) = S--(6 — 9b) for 

kink — > kink, and antikink — > antikink, as well as S-\ (9 — 9b) = S y{9 — 9b) for kink — > 

antikink, and vice versa. These S matrix elements are given by taking the massless limit 
of the results of H : 



c (q\ exp(A#) 

^ )= l + *exp(A0) ex P^)] 

s+ - {e) = iTi^W) eMiavm 



(4.1) 



Here exp[ia„] is the phase of the expression given in Eq.(3.5) of [JSj. For convenience we 
have defined 

A.I-!. 

V 

The boundary £ matrix is unitary: |£+_|_| 2 + \S-\ | 2 = 1. The S matrix is such that for 

particle with very large energy (UV limit), the scattering is diagonal. Diagonal scattering 
in the unfolded theory is totally off-diagonal scattering in the folded one, so matrix elements 
here are interchanged as compared to |§. 

Because we are no longer working in the plane- wave basis, it is also necessary to 
find the non-trivial S matrix for particles interacting in the bulk. Since all particles are 
massless, they must all have the same velocity vp (which is set to 1 in this paper). Thus 
it is not immediately obvious how to define an S matrix for two left movers. It is best 
interpreted as a matching condition on the Bethe wave function, as we will explain in 
the next subsection. Alternatively, one can think of acting on multiparticle states with 
creation and annihilation operators: the non-trivial S matrix means that these operators 
satisfy non-trivial commutation relations, the so-called Zamolodchikov-Faddeev algebra 
|24|1 . Since there is no mass scale in the bulk (only Tg at the impurity) , a two-particle bulk 



S matrix element can only depend on the ratio of the two momenta Pi/p2- In terms of 
rapidity, this is a function of 9\ —92- Thus we can now think of the impurity as a particle 
with rapidity 9b and a different S matrix than the bulk one. 

For general v, the left-left two-particle S matrix is given by the same formula as 
in the massive case. This S matrix is not diagonal: the initial state \K(9i)A(92)) can 
scatter to \A(9i)K(92)) because the kink K and antikink A have the same mass. For the 
kink/antikink scattering one has three amplitudes |24| 

a(9) = sin[A(7r + z0)]Z(0) 

b(6) = -sm(i\6)Z(0) (4.2) 
c(9) = sm(Xn)Z(9). 
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where the S matrix element a{6\— 62) describes the process \K(9i)K(9<2)) — > \K(6i)K(62)), 
as well as \A(9 1 )A(9 2 )) -> \A(6 1 )A(0 2 )), b describes KA -> KA, c describes the non- 
diagonal process KA —>■ AK, and there is a symmetry under interchange of kink to antikink 
(corresponding to <p° — > — <fi°). The function Z{9) is a normalization factor, which can be 
written as 

7{ n, 1 (• [°° *v ■ 2ex y sinh [( A - M \ 

Z{9) = — r—. —-t exp [ i — sin . 

sm|A(7r + %9)\ \ J_ 00 2y tv sum y cosh |AyJ / 



The breather-kink and breather-breather S matrices are well known [ 23 ] ; we do not write 
them down here. 

When 1/v is an integer, the bulk scattering is diagonal (c vanishes) and a = ±b. 
Therefore, the only allowed processes are of the form \j(0i))® |&(#2)) - * ^(^2))® 
Such a process is described by the S matrix element S^ lk {9i — 62). However, the impurity 
scattering is not diagonal, so charge transport is still possible. The bulk diagonal scattering 
makes the thermodynamic Bethe ansatz computation discussed in the next subsection 
much simpler at these values of v. 

4-3. The non-fermi distribution functions 

As with the S matrix, it is convenient to use rapidity instead of momentum or energy 
to define the densities of states and distribution functions. The number of allowed kink or 
antikink states per unit length with rapidity between 9 and 9 + d6 is given by n±(9, V)d6, 
while the number of states actually occupied by kinks or antikinks in this rapidity range 
is n+{9, V)f+(9, V)d9 and n_(9, V)/_(0, V)d9, respectively. Thus 

/oo 
d9 [n+(9,V) f + (9,V)-n-(9,V)f-(9,V)]. (4.3) 
-00 

Since at most one kink is allowed per level, we have < f± < 1. The functions rij and fj 
are defined for the breathers in the same manner. 

These thermodynamic functions rij(9) and f±(9, V) are different from the free-fermion 
functions when the odd-boson kink theory is an interacting Luttinger liquid [y 7^ 1/2), but 
we can derive them exactly. The idea is simple, and has become known as the thermody- 
namic Bethe ansatz |nj . It is basically the same as what is used in exact solutions of other 



impurity problems like the Kondo problem [20] . The main difference between the approach 



used here and the traditional approach is that in the latter, the Bethe ansatz appears as 
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a tool to diagonalize a bare hamiltonian, while here we work directly in the renormal- 
ized theory where factorized scattering is assumed, and the Bethe-ansatz equations (here 
relations between n and /) follow simply from a matching condition on the Bethe wave 
function. This matching condition gives one set of equations relating the functions nj and 
fj. Following the standard thermal approach to Bethe ansatz [ID], one writes the free 
energy as a functional of n and / and minimizes it. Together with the first equation from 
the matching condition, this second set of equations from the minimization requirement 
yields n and /. 

The simplest situation occurs when the scattering is completely diagonal, which occurs 
when 1/u is an integer. Here we can easily impose periodicity of the Bethe wave function 
since we know the S matrix, which encodes the change of the wave function when two 
particles (or a particle and the impurity) are interchanged. Thus we know the phase the 
wave function picks up when a particle is brought "around the world" through all the 
others. We denote generically the kink, antikink and breathers by the index j, so the 
two-particle S matrix elements are labeled by Sj k (9\ — 9<i). I n the following calculation 
of the densities n and /, we can ignore the phase shift due to scattering on the impurity 
because it affects densities only to order l/l, where I is the system size. This does not 
affect the densities in the thermodynamic limit I — > 00, and will therefore not show up 
in the computation of the conductance. (On the other hand, however, in the absence of 
the voltage these l/l effects on the densities can be computed and are responsible for the 
impurity free energy 0.) The requirement of a periodic boundary condition under this 
operation quantizes the momenta: 

Af 

e -im r . exp(e t )i/n "Q s r . r .(0. _ Q.} = 1 5 (4.4) 

where E = —P = rrij exp(6>) for a particle of type j. The periodicity condition ( fO| ) 
includes the interaction of any particle with all the others. By taking a derivative with 
respect to the rapidity of the logarithm of (4.4), one finds 

m p e r°° 

nj (e) = !!h^ + J2 de% k (o - e')n k {d')f k {e'), (4.5) 

11 u J-oo 



where $jfc(0) = -i(d/d9) In Sj^ lk (9) /2tc [0. For example, for v = 1/3 there is one 
breather (b) and p8| 

$ bb (9)= 2$ ++ (0) = 2$ + _(#): ' 



TTCOShtf 

^ <n\ ^ //n V2cosh6> 
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while the others follow from the symmetry + <-> — . Explicit expressions for $jfc(#) in the 
more general case where 1/u is a positive integer are given in Eq.(4.9) of 0, and in 
One defines an auxiliary pseudoenergy variable ej to parametrize fj via 



1 + e-^/ T e e i ' 



where the /x^ are the chemical potentials: /i + = — \i- = eV/2; /jj, = 0. By demanding that 
the free energy at temperature T (expressible in terms of fj and rij) be minimized, we find 
an equation for ej in terms of the (known) bulk S matrix elements: 



,0 roo 

^k/T-e k (0,V/T) 



m r°° 
ej (8, V /T) = ^--J2 d6 $> jk {6 - 6') ln[l + e 

k J -°° 



e-w> v ' r >\. (4.7) 



Solving this system of coupled integral equations for ej gives the functions fj. Except in 
special cases, this solution cannot be obtained in closed form, but it is easy to solve these 
equations numerically. We will not need the explicit because the equations ( |4.7| ) can 



be simplified by using an trigonometric identity described in [p9| , |30| for example. Denoting 
convolution by ★ 

f*g(6) = J de'f{6-d')g{9') 

one finds 

®jk(0) = N ^ K * CM*) + MW) ( 48 ) 
i 

where the kernel K{6) = \/{2-k cosh X9), and Nj k is the incidence matrix of the following 
diagram 

12 a A-2 / 

0--0 0---0--0 A-i 



i.e. Njk=l if the nodes j and k are connected, and otherwise (in particular Njj = 0). This 
identity follows simply by Fourier transforming the and using trigonometric identities. 
We show this explicitly for v = 1/3 (A = 2) in Appendix B. Thus the simplified form of 
( |yD is then 

€j(6, V/T) =K*J2 N 3k ln [l + e -^ k/T e ek ^' V / T) ] . (4.9) 

k 
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The dependences on the mass ratios seems to have disappeared from ( fO| ) , but they appear 
as an asymptotic condition: the original equations (|4.7|) indicates that the solution must 
satisfy 

rrij a 

€j — > — — e as — > oo. 
We emphasize for later use that comparing the relations (|4.7|) and ( |4.5|) gives 

n j (e,V) = ^d e e j (e,V). 
One effect of these equations is that the symmetry implies that 

n-(9,V) = n+(6,V) = n(9, V) 

and e_(0, V) = e + (6,V). 

The analysis at general values of 1/u is more complicated. Since the bulk S matrix 
is not diagonal, the phase picked up when bringing a particle around the world can be 
expressed only as an eigenvalue of a monodromy matrix, which itself must be diagonalized. 



This can be done at the price of introducing a further Bethe ansatz (see [|31| , |32[| for example) . 
The results are quite simple for vj (1—u) integer: the equation ([4.9D holds, with the addition 
of a term 5j\Me e /2 to the right-hand side. The general result is quite complicated, but 
the appropriate diagonalization has been done for all v J29| . We caution that only for \jv 



an odd integer does this one-component Luttinger liquid apply to the Hall effect; when v is 
a more complicated fraction, the edge states of the quantum Hall effect may be described 
by models with more than one boson 0]. 



5. The current and the conductance 

We have shown that the bosonic field theory with Hamiltonian (|3.6|) can be studied in 
terms of a particular set of quasiparticles and their scattering. We now compute an exact 
equation for the conductance using this basis. 

Without the backscattering, the left and right charges (or equivalently, the even and 
odd charges) are conserved individually. The backscattering allows processes where a 
charge carrier of the left-moving edge hops to the right-moving edge or vice versa. In the 
original basis, the current Is is the rate at which the charge of the left-moving edge is 
depleted. By symmetry, d t QL = —dtQn m each such hopping event, so Is = dt(|AQ) = 
d t yn^Q) 5 an d we see that in the even/odd basis, the tunneling corresponds to the violation 
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of odd charge conservation at the contact. In the S matrix language this happens when 

S.\ 7^ 0, so that a particle of positive odd charge (the kink) can scatter into one of 

negative charge (the antikink) at the contact. Neutral quasiparticles cannot transport 
charge and thus do not directly contribute to dfAQ. We emphasize that neither the bulk 
nor the boundary S matrix elements depend on the voltage; the voltage only affects the 
thermodynamic properties. 

To calculate the conductance, we start with a gas of quasiparticles with a chemical 
potential difference for kinks and antikinks corresponding to the voltage V. A positive 
voltage means that there are more kinks. When there are more kinks than antikinks, 
the backscattering will turn more kinks to antikinks than it turns antikinks to kinks. 
When a kink of momentum p is scattered into an antikink (the conservation laws require 
that it have the same momentum p) this changes AQ by —2. Since kink and antikink 
quasiparticles scatter off the point contact one-by-one, we may describe the rate at which 
this charge transport occurs in terms of two quantities: the probabilities of finding a kink or 

antikink of momentum p at the contact, and the transition probability | (9 — #b)| 2 - We 

therefore study the density of states n(9, V) = n+(9, V) = n-(9, V) and the distribution 
functions f±(9, V) in the thermodynamic limit (I — > oo) and in the presence of an applied 
voltage V. We can now compute the backscattering current from a rate (Boltzmann) 
equation. The number density of kinks of rapidity 9 which scatter into antikinks per 

unit time is given by | \ 2 n(9 : V)p^ , where p_| is the probability that the initial kink 

state of rapidity 9 is filled and the final antikink state is empty. In a fermi liquid, we 

would have simply p_\ = — /-]. However, in our interacting theory, correlations 

between the particles mean that this does not necessarily factorize in this manner. We 

can however write p_| = f + — , where /_| is the (unknown) probability that the 

kink and antikink state are both occupied. The rate at which antikinks scatter to kinks is 

likewise proportional to p |_ = /_ — /+-. Thus the charge Ql — Qr changes at a rate 

proportional to p_| — p |_ = /+ — /_, independent of the unknown factor. Using (|01 ) 

with I B = ed t Q°/V2 gives 

/oo 
d9 n(9)\S+-(9 - 9 B )\ 2 [f+(9, V) - /_((9, V)] . (5.1) 
-oo 

Without the backscattering, the current is Iq(V) = e(AQ)/2l, and it follows from ( fO| ) 
that Iq is given by the same expression without the — |^ | 2 . Thus the full current is 
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proportional to 1 — \S^ | 2 = IS^+I 2 . Using the definition of / and the result n(6>) 

^dge±(6) allows us to simplify the resulting expression. Notice that 

Defining 

e(0, V) = e+(0 - ln(M/2T), U) = e_(0 - ln(M/2T), U) 



and plugging into ( pTT|) gives our main result 

eT r°° 1 
I(T B ,V,T) = - I de i + eme _ 0B) de\n 



oo 

2/i 7.^ cosh 2 [A(0-ln(T B /T))] 



1 + e eV/2T e -e(9+ln(M/2T),V) 
1 _|_ e -eV"/2T e -e(e+ln(M/2T),"K) 
1 + e e^/2T e - e( e,y) 



1 _)_ e -eV/2T e -e(0,V) 



(5.2) 



where to get to the second line we integrated by parts and redefined 6 by a shift. Even 
though the breathers do not appear in ([5]^) , they interact with the kink and antikink and 
affect the calculation of e. 

The differential conductance is defined as G(Tb/T,V/T) = 8vI(Tb,V,T). In the 
V —>■ limit, the result can be written in a simple form. Using ( |4.7| ) and ( |4.5| ), it is easy 
to see that de/dV\v=o = 0, so 

e 2 A Z" 00 1 1 

G(T B /T,0) = ^ d9 —— = (5.3) 

1 ' ' 2h l + e^-o) cosh 2 [A(#-ln(T B /T))] 1 ' 

To check our result, we consider v = 1/2, where the conductance was previously 
derived exactly . As noted above, the odd-boson kinks are simply free fermions [p2|,|9|,p!9| , 
so they have the Fermi distribution function f±(9,V) = 1/[1 + exp((Me e =F e U)/2T)] 
implying that the function e(9) is simply t(Q) = e 6 . This of course follows from ( |4.7| ), 
because at v = 1/2 the bulk scattering is trivial and &jk = 0. Note again, as mentioned 
in Section 3.1 above, that we have used the kink/antikink description of the free fermion 
theory here, where all momenta have only one sign. Using the particle-hole transformation 
of Section 3.1, one finds immediately that the occupation number of kinks and antikinks 
are 

1 

f±(p) = ®(P) 1 + eipTV /2)/T i P > °> V = 2 

whereas the occupation number of unoccupied kink/antikink states are 

[i-/±](p) = e(p) 1 + e( _ p 1 ±W2)/T , P >o, v = \ 
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Now all momenta have one sign. Note that f±(p = 0,V = 0) = 1/2. These fermions (alias 
kinks) do scatter non-trivially off of the point contact, with S matrix given by ( f4.1| ). The 
resulting expression for G(T B , 0) obtained from (|5.3| ) is identical to the result in sect. VIII 
of H . Actually one can re-express the integral in terms of dilogarithm functions after some 
lengthy but straightforward manipulations to find 



I(T B ,V,u=±) 



e 2 V 
2h 



1 —I ruth I — h 

eV * \2 2tvT 



(\ , T B id 



+ 



4ttT J 



(5.4) 



where the digamma function ip(x) = F'(x)/V(x). At zero voltage in particular one finds 



e 

2h 



1 - 



-B 



2nT 



2 2tvT J 



(5.5) 



and at zero temperature 



I(V> 



e 2 V eT B eV 
arctan 



2h 



h 



2T, 



B 



(5.6) 



These v = 1/2 formulas also apply for v near 1/2 in the leading-logarithm approximation, 
after a V- and T-dependent renormalization of Tg B3 . 



Only at T = can these equations be solved in closed form for all T B ; we discuss this 
limit in the next subsection. However, we can study the solutions in certain limits. As 
T B /T — > 0, we can evaluate the conductance explicitly, The linear- response conductance 
( p73| ) becomes 

G(0, 0) = [/±(-oo, 0) - /±(oo, 0)]e 2 /h 

in this limit. One finds f±(oc,V) = obviously from ( [4.7| ) in this limit. More generally, 
one has 



eT 

J(Tb = 0,V,T) = -— In 



1 



:i-e 



iV/T 



)f+(-oo,V) 



(5.7) 



To find /(0, V), we use a well-known trick, given for example in [ i~l|j . The kernels in 
or ( [4.9| ) are appreciably different from zero only when 9' is near 9. Thus when we are 
interested in values of 9 near — oo, we can replace the value of e(9' , V) in the integral with 
e(— oo, V). We then can do the integral over the kernel explicitly. This then gives us the 
coupled difference equations 



N 



e j (-oo, V) = ln ( X + e-^ fe/T e £fe( 



) 
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Solving these explicitly gives /(— oo,0) = v and we indeed recover Go = ve 2 jh. For 
the non-equilibrium conductance, it is not even obvious that Go will not depend on V. 
However, one can check that 

,'eV(l-v)\ sinhz/eW2T 

f+{-oo,V) =exp 



IT J sinhel//2T 

and plugging this into ( |5"77| ) gives Go = ve 2 jh for all V. 

These non-perturbative equations also give the perturbative exponents. As Tb/T — > 
oo, the linear-response conductance G oc (T /T B ) 2{1 ~ v),v . Thus it goes to zero with the 
correct exponent, as in [[|. For Tb/T small, it can be argued that e(9 + in / (1 — u)) = e{9) 
p0(1 . Using this to write a power series for e and plugging into ( |5.2| ) gives G — Go oc 
(Tb/T) 2 ( 1-i/ ). Both are in agreement with ||. In fact, when v < 1/2 all of the coefficients 
g n {V/T) in the series G = Y^=o 92u(Tb /T) 2n ^ 1 ~^ can be computed using Jack polynomial 
technology Moreover, there is a non-perturbative functional relation [III relating the 
linear-response conductance to the free energy for all v < 1. 

We also note that there has been some confusion about the power of the exponent for 
T/Tb near zero. As discussed in ||, the leading (irrelevant) perturbation is the energy- 
momentum tensor, which is of scaling dimension two. This in fact means that the leading 
correction to the free energy is of order T 2 . (This power was also derived by scaling 
arguments in |[35|| .) However, this operator does not give any contribution to the DC 
conductance. Intuitively, this is because the energy-momentum tensor has no charge and 
should not affect charge transport. More precisely, one can check that when inserted into 
the Kubo formula, powers of the frequency u appear. When we take the DC uj — > limit 
this contribution vanishes. Thus the naive scaling arguments of [[35] do not apply to the 
DC linear-response conductance. However, outside of linear response at V > 0, one finds 
indeed that 

I(T B , V, T) w I(T B , V, 0) + T 2 I 2 (T B , V) + ... 



in agreement with p4]j35[ 



We can easily solve for the conductance numerically. To plot the complete function, 
one fixes a value of V/T and solves (|4.7| ) numerically for e to double-precision accuracy and 
inserts the result into ( |5.2| ). Evaluating the integral numerically for various values of T/Tb 
then gives I(Tb/T, V/T)/V as a function of T B for fixed V and T. To find the conductance 
it is easiest to just vary the voltage slightly in order to take the derivative numerically. A 
more precise way would be to use the fact that once one knows e numerically, the integral 
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equations for de/dV are linear and can be solved by inverting large matrices (of size the 
number of lattice sites used in the discretization of the integral) numerically. 

Several graphs of G for v = 1/3 are given in the figures. For V/T < 8, the plot is 
qualitatively the same as for V = 0: a flat region at G = e 2 /3h for T large, a transition 
region where the power series corrections cause it to fall off until it reaches its asymptotic 
form proportional to (T/Tb) 4 . However, at V/T ~ 8, a qualitatively new feature appears: 
G has a peak! It is not very sharp; the highest it gets for ^ = l/3isG~ .35e 2 //i. A 
variety of values of V/T are plotted in fig. 1 and fig. 2. In the first we plot it versus Tb/V 
in order to make the approach to the T = limit to be discussed in sect. 6 clear. In the 
second we plot it versus Tb/T; we see that there is a substantial broadening of the curve 
as the voltage is increased. This should provide a prominent signal in the experiments. 

This peak and the values of V/T for which it occurs can be understood theoretically. 
A peak occurs for voltages large enough so that g 2 (V/T) changes sign. This must happen 
because at zero temperature the current can be expanded for large voltage in the form 

rp \ 2(1-1/) 

J-B 



l(y)<xV(v + C{-f) +...) 

Because the backscattering must make the current decrease we have C < 0. Taking the 
derivative with respect to V we see that for v < 1/2 the conductance must increase for 
small enough V/Tb- We can even find the value V* where g 2 (V*/T) = analytically, 
because g 2 can be calculated in closed form [ |T2"|| . This is done by first calculating the 
coefficient Z 2 of the partition function in the case where the interaction (in Euclidean 
time) is cos [V20(O, r) + 2irprT\, which is @ 



,_ , . sin7r^ r(l — 2u) 

Z2KV) - 



sin7r(z/ + — v + p)T(l — v — p) 

The finite voltage case is obtained by analytically continuing 2rrp to iveV/T |J. Then we 
use the relation g 2 (V/T) oc RedpZ 2 (p) (which can be shown by using explicit perturbation 
theory in the impurity interaction) to give 



g 2 (V/T) cxRe 



iueV\ , ,. iveV \ , iveV ' 
{^{y + ——) - V(l - v + ——)} 



(5.8) 



2nT J 2nT ' v 2nT 

Solving for g 2 (V*/T) = gives eV* /T = 7.18868564374998... for u = 1/3 and eV* /T = 
6.653022289582846... for v = 1/4. 

The perturbative results of |T^] give a completely independent check on the TBA 
results. One can fit the TBA results to a power series numerically to determine the 
perturbative coefficients to compare; one indeed finds for example that g 2 (V* /T) = in 
the TBA results. One can also check that the relation Tg = CA^ 1 v is independent of V, 
showing that the boundary S matrix is indeed independent of voltage. 
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6. Explicit Solution at T=0 

Some remarkable simplifications take place in the T = limit. Although the densities 
of states are still non-trivial, the distribution functions / become step functions. As a re- 
sult, the TBA equations become linear and can be solved using the Wiener- Hopf technique. 
We find explicit series expressions for I(V/Tb) and G(V/Tb) in this limit. Moreover, this 
leads to an exact duality between large V/Tb and small V/Tg. In the Hall devices this 
corresponds to a duality between Laughlin-quasiparticle tunneling and electron tunneling. 

At T = and V = the ground state of the theory is just the vacuum with neither 
kinks nor breathers; these particles are in fact defined as excitations above this vacuum. 
When V is turned on, this ground state becomes unstable. For V > 0, kinks of charge e 
start filling the vacuum, since they are energetically favorable for small enough momentum 
(large negative rapidity). The new ground state is made of kinks occupying the range 
9 E (-oo, A]; in other words, f + (9, V) = 1 for 9 < A and f+(9, V) = for 9 > A. The 
surface of the sea is approximately A m ln(eV/M), but computing A exactly requires some 
technology because the kink interaction affects the filling of the sea. There are no antikinks 
nor breathers in the sea at T = 0, so their densities do not appear in this analysis. For ease 
of notation we define p(9) = n+(9, V)f+(9, V)\t=o- When T = the periodicity relation 
( fi~5l) reduces to the following coupling between kink rapidities: 



where p = in [A, oo), and $ = doS^+^O) /2ix follows from the kink-kink bulk S matrix 



where in the above example g(9) = ffe e for 9 G (—oo, A], g(9) = otherwise. By taking 
Fourier transforms, and since p vanishes outside (— oo, A] we have 




(6.1) 



We consider now the general equation 




(6.2) 




(6.3) 



where we defined Fourier transforms by 
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For any function h we also introduce 



h- (u) = h(u})e 



-iujA 



This subscript should not be confused with the antikink index used in other sections of 
this paper. 

The integral equation ( |6.2p for p can now be written in the form 



p{l-*)-g = e vaA X + [u) i (6.4) 

where X+(u) is analytic in the upper half plane. Since p vanishes in [A, oo) we see that p- 
is analytic in the lower half plane, and similarly for gr_. To implement the Wiener- Hopf 
technique (see the appendix in |36] for a very clear exposition) , we need to factorize 1 — $ 



into the form 

l-$ = dnh[(A + lW2A] (65) 
2cosh[7ru;/2]sinh[7ru;/2A] G + G_ ' v ' ; 

where G_|_ (G_ ) is analytic in the upper (resp. lower) half plane and G_|_ (G_ ) vanishes only 
in the lower (resp. upper) half plane. One finds 

r /9W> i ri rH(i + AV/2A) „ A 

G + = V2tt(A + 1) r( _ iw/2A)r(1/2 _ iw/2) e > (6-6) 

G_(w) = G+(-u;), 

A = ilnA-i^ln(l + A). 

The phase ensures analyticity of G_|_ at zoo. Having set up the equation in this form, we 
can write down the solution for p: 

^- = [g-G+U (6.7) 



where 



In our case 



so we find 



2«7r j-oo u/ — u; + iO 



2/i + 1 



HK J 2ih u-i K 
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This relation was obtained also in [[57 . 

We have not yet determined the dependence of the cut-off A on the physical potential 
V. This can be done in two ways. The first way consists in minimizing the energy of the 



system. Following [35] we introduce a function e(9) which satisfies 



eV _ M 6 

~T~Y e 



e(9) 



$(6 - 6')e(6')d6' , 



(6.9) 



and the cut-off follows now from the condition e(A) = 0, which reads in terms of fourier 
transforms 

lim uje-(iuj) = 0. 

UJ — >oo 

We can easily solve for e. This is the same equation as the one for p except for the 
substitution 

9(9) = ^ - -e e , ~g(u)=ireVd(u)-- e ——-. 
2 2 2 iuj + 1 

We find therefore 



, , MG-(u)G + (i) A eV G-(u;)G + (0) 
e-(u) = - — w e A + 



2i 



uj — i 



2i 



(6.10) 



It follows that 



eA = eV G+(0) 



(6.11) 



M G+(i) 

The second method of determining A is to study the energy in the absence of backscat- 
tering, which is 



E 



d9p(9) 



M g _ eV 
T e ~ ~2~ 



- f «o) 



M2 -G + (i)G-(-i)e 2A - ^G_(0)G + (i)Ve 



A 



8h 



Ah 



This energy results from the massless description obtained by considering the free field 
as the A — > limit of the sine-Gordon model (|3.6|) . We can recover this energy directly 
from the Hamiltonian in this limit. Coupling to a potential amounts to considering the 
Hamiltonian (in the original left-right basis) 



h f l eV f l 

Hn = - — / dx [(d x (p R ) 2 + (d x (p L ) 2 ] - — I dx [d x (p L + d x (p R ] . 



4tvu 



(6.12) 
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Redefining the fields 4>l and (p B by a shift removes the linear term, and gives the energy 
per unit length associated with the potential V. Equating the two energies gives A in 
agreement with formula (|6.11| ). 

Since we know the density explicitly from ( p.8|) and ( |6.11| ), an explicit form for the 
current follows from ( |5.1|) , which at T = reads 



(6.13) 



I(V,T B ) = e / p(8)\S ++ (e-e B )\ 2 d6 

J — oo 

e 2 v r° . 1 



de f{9) 



2h v 1 + exp[-2A(# + \n(eV/T' B ) - A)] ' 

where 

p s G-MG + (0) 
1 + iu 

and the boundary temperature is redefined: 

rpl m ^ —A G+W rp g V^t A + 1 ) T( \ + ± ) 

= G^0)- TB T(X) ' 

We can therefore extract the following results from (JfTTtJ). At very large potential, 
the S matrix element goes to one and Iq = GqV, with 

e 2 r° e 2 u 

G " = Ykij (e)M = ir- (fU4) 

For small voltage, we can expand the S matrix element in powers of V/T' B to get 



2t/ /-0 00 / rr \ ' 

™ =4 L F(e)de ^t )n+1 



n=l 

where 



, N (6.15) 

fv^ T ( <V 



2n 



2tt 



f g-(^)g+(0) , = f +1 G_[-2mA]G + (0) 
7.^ (1 + iw)(2nA - iu) 1 J 1 + 2nA 



This expansion is valid only for eVe A /T' B < 1. By using (|6.6|) , we find our final result 
for the low- voltage expansion 
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We prove in the appendix that this strong-barrier expansion can be transformed into 
the following weak-barrier expansion 



_ e 2 Vv 
h 



oo , — -p/ \ / -it\ 2n(u— 1) 



1-2JC- 1 ) 

n=l 



n+1 u^T(nu) ( eV 



2T(n)m+n(u-l)) \T' B 



T^a> 1 - ( 6 - 17 ) 



Up to a constant and a shift, this is the same expression as for strong backscattering, with 
the replacement v — > We therefore have the result 

1(7%, V,u) = ^- v 2 I(T' B , V, -). (6.18) 
h v 

(In this equation, T' B is treated as a constant — it is the same on both sides, although 
we have seen T' b /Tb depends on v.) This proves that for T = at least, the weak- 
barrier (small Tg) and strong-barrier (large Tg) are completely dual to each other. In the 
Hall device, this duality is a concrete proof of a long-suspected relation between electrons 
and Laughlin quasiparticles. This arises from the fact that for weak backscattering, the 
operator for the tunneling of Laughlin quasiparticles is relevant while the electron tunneling 
operator is always irrelevant (|||. Laughlin quasiparticles are allowed to tunnel because 
they are tunneling through the bulk of the sample. However, in the strong backscattering 
limit, only electrons can tunnel because they are not tunneling through the Hall fluid. Thus 
the least irrelevant operator arises from electron tunneling. As seen in ||, the associated 
exponents are indeed related by duality; we find here that the entire expansions around 
these two limits are related. 

This duality was somehow known in the literature [[39]], but its status is not totally 
clear to us. It is usually considered as only approximate since it relies on an instanton 
approximation in the large-barrier limit 0. However, one can prove an exact duality |10 



(in the context of dissipative quantum mechanics, which is equivalent) between the cosine 
problem in the UV and the "tight-binding" problem in the zero-temperature limit, and 
then in [^IJ the tight-binding problem is mapped back onto the cosine problem, providing 
a sort of proof. The main source of difficulty is that the zero-temperature action must 
be handled with great care. Perturbation theory around the zero-temperature fixed point 
is ill-behaved and depends on an infinite number of counterterms (see the discussion at 



the end of |y] in the case of the flow from tricritical Ising to Ising model), so identifying 
the leading term in the approach to this fixed point is not sufficient. This is equivalent 
to saying that what one calls the strong-barrier problem must actually be defined with 
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great care. The strong-barrier problem which is at the end of our renormalization-group 
trajectory follows formally from dimensional continuation of the integrals for the weak 
barrier problem. It is not in any case a generic strong-barrier problem. For instance 
regularizing the integrals with a short-distance cut-off would give very different results, 
with a non-monotonic conductance []40j . Also, note that the duality does not apply to the 
impurity free energy, for example as explained in for small T the leading contribution 
is proportional to T 2 for all v. 

Interestingly ( |6.17| ) can be compared with the perturbative expansion in || 



2 



(6.19) 



thus providing the relation between the parameter Ai in the action and the TBA parameter 

AiK"" = ^T(v) (vT'z) 1 -. (6.20) 

where k is the non- universal (dimensionful) cutoff which appears in the boson two-point 
function, as defined in |12]| . 

Finally, setting v = ^ (i.e. A = 1) in these formulas gives G+(m) = 1 so I n = 2 n+i 
and T' B = 4Tb- Hence 

e 2 V ^ (-1)" +1 feV\ 2n e 2 V eT B eV 

arctan 



y ' 2h ^ 2n + l \T B J 

n=l x ' 

in agreement with the T —>■ limit of eq (|5.4|) 



2h h 2T, 



B 
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Appendix A. Duality 

In this appendix we provide details of the derivation of duality between strong- and 
weak-tunneling regimes. Defining 

eV 

we see that there are two different regions to consider in the computation of the integral 



1 + 66-2(^-1)6 



depending on whether be~ 2 ^~ l ^ e is greater or smaller than one. In each region, one can 
expand the integrand appropriately, and integrate term by term. For V/Tb small, this 
piece is always smaller than one, and the expansion ( |6.15|) follows immediately. For V/Tb 
large, one divides the integral in (]6.13|) into two pieces to get 



2y 

HV) =^PP 



du G_(u)G+(0) 



2n 



1 + ioo 



Em: 



N+l 



<2N 



.N=0 



+6 -iw/(3(i-l)) 



N=0 



+ 



iuj + 2(±-l)N 
(-1)" 



_iu + 2(±-l)N -iw + 2(£-l)(JV + l)_ 



(A.l) 



The poles of the rational functions except the one for ui = are all cancelled out because 
G-(2N(^ — = for N > 0. For the first sum we close the contour in the lower half 
plane. Since G_ is regular in this regime, the only pole is at uj = 0. For the second sum 
we close the contour in the upper half plane where there are poles of G-{uj) at 



u = 2in(u — 1), n > 0. 
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as well as those of the rational function at u = 0. Collecting all terms gives 



I(V) 



e 2 V 



2h 



\2u + ^{y - 1) 7^ 

I l + 2n(i/-l ; 



T(-nu)T(l/2 + n(u- 1)) 



1 



, x 2n(l-!/) 



B 



2(i-l) Ve^ 
Using the identity 



7V=0 



-1 



+ 



(-1) 



A' 



- H. 

2-^i N - 



+ 



N-nu (N + l) + nu 



(A.2) 



7T 



7V=0 



N — nv N + l + nu 



smnnv 



together with standard gamma- function identities gives ( |6.17|) . 



Appendix B. Incidence Matrix 

In this appendix we give some details involving the manipulations of the "incidence 
matrix" of section 4, for the case v = 1/3. We start by writing 





(®++ 


















$-6 





Using the explicit expressions given in ( |4.6|) , one finds that this may be written in the form 

*(0) = $ ++ (6>)N 2 + $ +b (0)N 

where 

N = 




It is convenient to fourier transform the rapidity dependence, in order to turn convolutions 
into multiplications. We find 

1 



/oo 
dde ike $ ++ (d) 
-oo 

/oo 
d6e lke <S> +b {6) 
-oo 



2cosh(fc7r/2) 



cosh(/c7r/4) 



cosh(/c7r/2) 

Making use of the relation N 3 = 2N, we find in fourier space 

N ■ *(fc) = $ +b (k)N 2 + 2$ ++ (fc)N 
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If we define the multiplicative factor 

1 



K{k) = $ ++ (k)/$ +b (k) 



2cosh(/c7r/4) 
we find 

KN # = $ + (2K&++ - $ +6 )N 

The expression in parenthesis equals (2KQ ++ — $+&) = 1/(2 cosh(/c7r/4)) = K(k). In 
conclusion we have found 

KN-$ = * + M (B.l) 

Fourier-transforming back to rapidity space, we find the relation (|4.8|) of section 4. Note 
that for v = 1/3, 

J-oo 27r 



7rcosh(20) 



as indicated below 
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Figure Captions 



Fig. 1 . The non-equilibrium conductance as a function of \og(Ts /V), for various values 
of V/T. Notice that the curve develops a peak when V/T > 7.18868, and that the T -> 
limit is smooth. 

Fig. 2. The non-equilibrium conductance as a function of Tb/T for various values of V/T. 
The curve broadens substantially as V/T is increased, even before developing a peak. 
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